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1.  OBJECTIVES 

The  objective  of  the  research  is  to  develop  efficient  computational  tools  for  non-linear  analysis 
of  smart  composites  structural  systems  with  embedded  smart  (e.g.,  piezoelectric  or 
magnetostrictive  particle)  layers.  The  main  thrust  of  the  contract  is  to  develop  analytical  as  well 
as  computational  (i.e.,  finite  element  analysis)  procedures  for  the  stress  and  vibration  response 
of  laminated  composite  plate  structures  with  smart-material  layers.  The  first  and  second 
phases  of  the  research  was  concerned  with  (1)  the  theoretical  formulation  and  analytical 
solutions  for  vibration  suppression  of  laminated  plates  and  (2)  the  linear  finite  element  analysis 
of  laminated  composite  plates  with  smart-material  layers,  respectively.  This  report  summarizes 
the  work  carried  out  in  the  third  and  last  phase  of  the  contract;  in  particular,  it  summarizes 
the  linear  transient  analysis  using  the  FEM  and  non-linear  finite  element  formulation  based  on 
the  third-order  shear  deformation  plate  theory. 


2.  INTRODUCTION 

The  study  of  smart  materials  and  structures  has  received  considerable  attentions  in  recent 
years.  The  advantage  of  incorporating  these  special  types  of  materials  into  the  structure  is  that 
the  sensing  and  actuating  mechanism  becomes  part  of  the  structure  so  that  one  can  monitor 
the  structural  integrity/ health  of  the  structure.  There  are  a  number  of  materials  that  have  the 
capability  to  be  used  as  a  sensor  or  an  actuator  or  both.  Piezoelectric  materials, 
magnetostrictive  materials,  electrostrictive  materials,  shape  memory  alloys,  and 
electrorheological  fluids  provide  examples  of  such  materials.  Among  these,  piezoelectric  and 
magnetostrictive  materials  have  the  capability  to  serve  as  both  sensors  and  actuators. 
Piezoelectricity  [1]  is  a  phenomenon  in  which  some  materials  develop  polarization  upon 
application  of  strains.  Examples  of  piezoelectric  materials  are  Rochelle  salt,  quartz,  and  Lead 
Zirconate  Titanate  or  PZT  (Pb  (Zr,Ti)  O3).  Piezoelectric  materials  exhibit  a  linear  relationship 
between  the  electric  field  and  strains  for  low  field  values  (up  to  100  V/mm);  and  it  exhibits 
nonlinear  behavior  for  large  fields,  and  the  material  exhibits  hysteresis  [2].  Furthermore, 
piezoelectric  materials  show  dielectric  aging  and  hence  lack  reproducibility  of  strains,  i.e.,  a 
drift  from  zero  state  of  strain  is  observed  under  cyclic  electric  field  conditions.  Terfenol-D,  a 
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magnetostrictive  material  [3],  has  the  characteristics  of  being  able  to  produce  strains  up  to 
2500pm  and  energy  density  as  high  as  25000  J /m3  in  response  to  a  magnetic  field. 


Vibration  and  shape  control  of  flexible  structures  is  achieved  with  the  help  of  actuators  and  a 
feedback  control  law.  Many  modern  techniques  have  been  developed  in  recent  years  to  meet 
the  challenge  of  designing  controllers  that  suit  the  function  under  some  required  conditions. 
There  have  been  a  number  of  studies  on  vibration  control  of  flexible  structures  using  smart 
materials.  Anjanappa  and  Bi  [4,5]  investigated  the  feasibility  of  using  embedded 
magnetostrictive  mini  actuators  for  smart  structure  applications,  such  as  vibration 
suppression  of  beams.  A  self-sensing  magnetostrictive  actuator  design  based  on  a  linear  model 
of  magnetostrictive  transduction  for  Terfenol-D  was  developed  and  analyzed  by  Pratt  and 
Flatau  [6].  Eda  et  al  [7]  and  Krishna  Murty  et  al  [8,9]  proposed  magnetostrictive  actuators  that 
take  advantage  of  the  ease  with  which  the  actuators  can  be  embedded  and  remote  excitation 
capability  of  magnetostrictive  particles  as  actuators  for  smart  structures.  Reddy  and  Barbosa 
[10]  presented  a  general  formulation  and  analytical  solution  for  simply  supported  boundary 
conditions  of  laminated  composite  beams  with  embedded  magnetostrictive  layers.  All  these 
studies  were  based  on  free  vibration  analysis  of  the  structures  with  velocity  feedback  control. 


In  the  present  study,  control  of  the  transient  response  of  laminated  composite  plates  with 
integrated  smart  material  layers,  used  as  sensors  and  actuators,  is  studied  using  a  unified 
plate  theory  that  includes  the  classical,  first-order  and  third-order  plate  theories  as  special 
cases.  A  simple  negative  velocity  feedback  control  is  used  to  actively  control  the  dynamic 
response  of  the  structure  through  a  close-loop  control.  A  displacement  finite  element  model  of 
the  equations  of  motion  governing  the  unified  theory  is  developed.  The  effects  of  material 
properties,  lamination  scheme,  and  placement  of  the  smart  material  layer  on  deflection 
suppression  are  investigated. 

3.  THEORETICAL  FORMULATION 

3.1  Introduction 

The  simplest  equivalent  single-layer  theories  are  the  classical  laminate  plate  theory  (CLPT)  and 
the  first-order  shear  deformation  theory  (FSDT).  These  theories  describe  the  kinematics  of 
most  laminated  plates  adequately.  The  third-order  shear  deformation  theory  (TSDT)  represents 
the  plate  kinematics  better  and  yields  better  inter-laminar  stress  distributions.  Quadratic 
variations  of  the  transverse  shear  strains  and  stresses  through  the  layer  avoid  the  need  for 
shear  correction  coefficients  as  required  in  the  first-order  theory. 

3.2  Displacement  field 

The  displacement  field  for  the  third-order  shear  deformation  theory  (TSDT)  can  be  expressed  in 
the  form  (see  Reddy  [12]) 


3  Swn 

u(x,y,z,t)  =  UQ(x,y,t)  +  z^x(x,y,t)-c\z  (tyx  +— — ) 

ox 

3  dwQ 

v(x,  y,  z,  t )  =  vo  (x,  y,t)  +  z$y  (x,  y,  t)-c\z  (<[> y  +  — ) 


w(x,y,z,t)  =  wq  (x,y,t) 


(1) 

(2) 

(3) 


where  (w0,v0,w0)  and  (<t>x,<l>y)  have  the  same  physical  meaning  as  in  the  first-order  theory; 
they  denote  the  displacements  and  rotations  of  a  transverse  normal  on  the  plane  z  =  0, 
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respectively.  Then  the  displacement  field  of  FSDT  is  obtained  by  setting  cl  =  0,  and  for  the 
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Reddy  third-order  theory  it  is  equal  to  cj  =  4/3/?  . 

3.3  Equations  of  motion 

The  equations  of  motion  of  the  third-order  shear  deformation  theory  are  derived  using  the 
dynamic  version  of  the  principle  of  virtual  displacements.  They  are  (see  Reddy  [12]) 
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and  (Nxx,Nyy,Nxy)  denote  the  total  in-plane  force  resultants,  [M  XX,M  yy ,MXy )  the  moment 
resultants,  and  (PXX’Ryy’Rxy)  and  ( Rx,Ry)  denote  the  higher-order  stress  resultants. 
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The  force  and  moment  resultants  are  related  to  the  strains  by 
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where  the  stiffnesses  A- ,  Z)  ,  and  Z.  are  defined  for  i,  i  =  1,  2,  6  as  well  as  i,  i  =  4,  5.  The 
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stiffnesses  B  ,  E.. ,  and  H  are  defined  only  for  i,  j  =  1,  2,  6.  The  coefficients  of  A.. ,  B ,  Dr  , 

ij  7  ij  7  ij  J  J  J  V  J 

Etj  ,  ,  and  Fly  are  given  in  terms  of  the  layer  stiffnesses  Q  and  layer  coordinates  zk+l  and 


The  stress  resultants  associated  with  the  magnetostrictive  materials,  {NM},  {Mm},  and 
{ P  M  }  are  defined  by 
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where  e„  is  the  transformed  moduli  of  the  actuating/ sensing  material,  which  in  the  present 

study  is  taken  to  be  a  magnetostrictive  material  ,  and  H  _  is  the  electric  field  intensity  which 
should  be  excluded  in  the  constitutive  relations  for  the  structural  part  of  the  composite 
structures. 

3.4  Velocity  feedback  control 

Considering  velocity  proportional  closed-loop  feedback  control,  the  magnetic  field  intensity  H 
can  be  expressed  in  terms  of  coil  constant  kc  and  coil  current  l{x,y,t )  as 


H(x,y,t)  =  kc  I(x,y,t)  (14) 

and 
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b :  +  4  r: 


I(x,y,t)  =  c(t)^~ 
ot 


where  bc  is  the  coil  width  ,  rc  is  coil  radius,  nc  is  number  of  turns  in  the  coil,  and  c(t)  is  the 
control  gain  which  is  assumed  to  be  a  constant  in  this  study. 

4.  FINITE  ELEMENT  MODEL  (LINEAR) 

4.2  Virtual  work  statements 

The  weak  forms  of  the  equilibrium  equations  are 
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where  V ,,  is  defined  as 
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The  primary  variables  of  the  third-order  theory  are  uu,us,wq, - ,(])„, ())5,  where  (un,us) 

dn 

denotes  in-plane  normal  and  tangential  displacements,  and  (<|)n,<t)s)  are  the  rotations  of  a 
transverse  line  about  the  in-plane  normal  and  tangent.  Lagrange  interpolation  functions  of 
{lln,u s,(j)n,(j)s)  and  Hermite  interpolation  function  for  VV0  are  used  for  the  formulation  of  the 
displacement  finite  element  model.  A  conforming  element  that  has  eight  degrees  of  freedom 

(«0  >  V0  >  W0  ■>  W0,x  >  W0 ,y  ’  W0,xy  >  </>x  >  </>y  )  is  USed  in  this  Study- 


4.3  Semidiscrete  finite  element  model 


The  generalized  displacements  are  approximated  over  an  element  Q1  by  the  expressions 
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where  p ?  denotes  the  Lagrange  interpolation  functions  and  cp?  are  the  Hermite  interpolation 
functions.  Here  we  chose  the  same  approximation  for  the  in-plane  displacements  (u0,v0)  and 
rotations  (<|)x,<|>j, ),  although  one  could  use  different  approximations  for  these  two  pairs.  In  the 
case  of  the  conforming  element,  the  four  nodal  values  associated  with  vv0  are 
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The  finite  element  model  is  of  the  compact  form 


6 


Z  Z  (KfA 7  +  C?  AV  +  A^)-  Fta  =  0  ,  I  =  1,  2, 

P= i  i=i 
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where  dc{t )  =  e mdmk cc(t)  and  em  is  modulus  of  the  magnetostrictive  layer,  dm  the  magneto¬ 
mechanical  coupling  coefficient,  and  kr  and  c(l  )  are  defined  in  the  previous  section.  For  the 
details  of  the  stiffness  and  mass  coefficients,  see  Reddy  [12]. 

This  completes  the  general  linear  finite  element  model  development  of  the  third-order  shear 
deformation  plate  theory. 

4.4  Transient  analysis 

The  equations  of  motion  can  be  solved  exactly  using  analytical  methods,  but  those  are 
algebraically  complicated  and  require  the  determination  of  eignevalues  and  eigenfunctions,  as 
in  the  state-space  approach.  Newmark  method  that  takes  advantage  of  the  static  solution  form 
for  spatial  variation  and  uses  a  numerical  method  to  solve  the  resulting  differential  equations 
in  time  was  used  to  determine  the  transient  response  of  composite  laminates  in  this  study. 

The  second-order  differential  equation  (27)  can  be  expressed  in  matrix  form  as 

[Me]  {Ue}  +  [Ce]  {id}  +  [Ke]  {ue}  =  {Fe}  (33) 

Using  the  Newmark’s  scheme,  Eq.  (33)  can  be  reduced  to  the  form 

[K]s+l{A}s+l={F}s;s+1  (34) 
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where 


[KU=[K]s+l+a3[M]s+l+a6[C]s+l 

{F}s^={F}s+l+[M]s+l{A}s+[C]s+AB}s 

{ As  =aMs  +aMs  +a5 R  ,  {^}s  =  a6  {u}s  +  a7  {ii}s  +  as  {ii}s 

a,  =  (1  -  a)At ,  a2  -  aAt ,  a2  =  — — — ,  a4  =  a,  At ,  a5  =  — — — 

r{ AO  '  r 


2a  ,  .  . 

a  7  = - 1 ,  ag  =  A? 


A 


(35) 

(36) 

(37) 


(38) 


and  a  and  y  are  parameters  that  determine  the  stability  and  accuracy  of  the  scheme. 

In  this  numerical  method,  the  time  derivatives  are  approximated  using  difference 
approximations,  and  therefore  solution  is  obtained  only  for  discrete  times  and  not  as  a 
continuous  function  of  time. 

5.  NUMERICAL  RESULTS  AND  DISCUSSION 
5.1  Preliminary  comments 

Numerical  studies  are  carried  out  to  analyze  the  deflection  suppression  characteristics  for 
different  lamination  schemes  and  boundary  conditions  using  the  developed  TSDT  finite  element 
models.  The  baseline  of  the  simulations  is  the  simply  supported  square  laminate  (a/b  =  1, 
a/h  =  10;  see  the  Figure  1)  under  sinusoidal  distributions  of  the  initial  velocity  field 


dw 

~dt 


,  rv\  •  TOC  .  TCV 

(x,  y,t  =  0)  =  sin  —  sin  — 
a  b 


(39) 


The  time  step  selected  in  the  present  study  is  t  =  0.0005sec.  The  notation  for  lamination 
scheme  (0l,d2,03,d4,m)s  means  that  there  are  10  layers  symmetrically  placed  about  the 

midplane  with  the  fiber  orientations  being  (0l,02,03,04,m,m,04,02,02 ,0{),  where  in  stands  for 
the  magnetostrictive  layer  and  subscript  S  stands  for  symmetric.  The  material  properties  of 
Terfenol-D  and  the  elastic  composite  materials  are  listed  in  Table  1. 

In  finite  element  analysis,  solution  symmetries  should  be  taken  advantage  of  by  identifying  the 
computational  domain  to  reduce  computational  effort.  For  a  laminated  composite  plate  with  all 
edges  simply-supported  or  clamped,  a  quadrant  of  the  plate  may  be  used  as  the  computational 
domain.  Figure  2  shows  two  types  of  simply  supported  boundary  conditions  for  the  third-order 
shear  deformation  theory.  Figure  3  shows  the  effects  of  the  finite  element  results  for  the 
different  laminations  and  boundary  conditions.  Quarter  plate  models  with  proper  boundary 
conditions  can  be  used  in  the  antisymmetric  laminates  with  simply  supported  boundary 
condition,  but  not  with  the  clamped  cases.  For  symmetric  laminates,  the  simply  supported 
cross-ply  laminates  can  be  modelled  as  a  quarter  plate.  The  boundary  conditions  along  a  line  of 
symmetry  should  be  correctly  identified  and  imposed  in  the  finite  element  model.  When  one  is 
not  sure  of  the  solution  symmetry,  it  is  advised  that  the  whole  plate  be  modelled. 


8 


Table  1  Material  properties  of  magnetostrictive  and  elastic  composite  materials. 


Material 

Ei 

[GPa] 

e2 

[GPa] 

Gi3 

[GPa] 

G23 

[GPa] 

Gi2 

[GPa] 

V'n 

P 

[Kg  nP3] 

7 

O 

Terfenol-D 

26.5 

26.5 

13.25 

13.25 

13.25 

0 

9250 

1.67 

CFRP 

138.6 

8.27 

4.96 

4.12 

4.96 

0.2 

6 

1824 

- 

Gr-Ep(AS) 

137.9 

8.96 

7.20 

6.21 

7.20 

0.3 

0 

1450 

- 

Gl-Ep 

53.78 

17.93 

8.96 

3.45 

8.96 

0.2 

5 

1900 

- 

Br-Ep 

206.9 

20.69 

6.9 

4.14 

6.9 

0.3 

0 

1950 

- 

1234  5678  9 

Equivalent  8x8  Full  Plate  Model 


Fig.  1.  Laminated  composite  plate  configurations  and  finite  element  meshes. 
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SS-l  conditions;  v0=  wQ  =  <j>  =  0 
Nxx=  Mxx=y° 

SS-2  conditions;  u0  =  w0  =  (J)y  =  0 

Nxy=Mxx=0 


Simply  supported 
boundary  conditions 


SS-l  conditions;  u0  =  w0  =  <J) .  =  0 
Nyy=M„=0 

SS-2  conditions;  vo  =  wo  ~  =  ^ 

N  =  M=  0 


x 


Fig.  2. 
plates. 


Two  types  of  simply  supported  boundary  conditions,  SS-land  SS-2,  for  laminated 


Time  (sec) 


Time  (sec) 
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(c) 


Time  (sec) 


/  i  v  Time  (sec) 

(d) 

Fig.  3.  Effects  of  using  full  and  quarter  plate  models  in  finite  element  modelling  of  CFRP 
composite  plates;  cross-ply  (m,90,0,90,0)SOranti-s,  angle-ply  (m,30,-30,30-30)s  or£mti-s,  and  general 
angle-ply  (m,45,-45,90,0)s  or  anti-s  ;(a)  Simply  supported  antisymmetric  laminates,  (b)  Clamped 
antisymmetric  laminates,  (c)  Simply  supported  symmetric  laminates,  (d)  Clamped  symmetric 
laminates. 

5.2  Simply  supported  laminates 

To  compare  with  the  analytical  results,  the  SS-1  boundary  conditions  and  quarter  plate  model 
have  been  used  for  symmetric  cross-ply  laminates.  The  deflections  predicted  from  the 
analytical  (eigenvalue  analysis)  and  transient  finite  element  analysis  are  within  the  reasonable 
agreement,  as  shown  in  Figure  4.  Figure  5  shows  the  central  displacements  using  the  different 
plate  theories  (CLPT,  FSDT,  and  TSDT)  for  two  different  lamination  schemes.  It  is  observed 
that  CLPT  shows  higher  deflection  suppression  capacity  in  both  cases.  This  is  expected 
because  the  CLPT  renders  plate  stiffer  compared  to  the  other  theories  (also  see  Table  2). 

After  studying  the  influence  of  lamina  material  properties  on  the  amplitude  of  vibration  and 
vibration  suppression  times,  it  is  observed  that  materials  having  the  almost  same  EjE7  ratio 
have  similar  vibration  suppression  characteristics.  Figure  6  shows  the  deflection  damping 
characteristics  for  the  different  laminate  materials. 
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0.010 


0.005 


g  0.000 

o 

ft 

CO 

5  -0.005 


-0.010 


Analytic  solution  - FEA  solution  (dt=0.0005) 

-  FEA  solution  (dt=0.0001)  - Uncontrolled  motion 

A  A  A  A  A 

’  /  \  '  V  /  \  /  \ 

I  \ 


0.00 


0.05 


0.10 

Time  (sec) 


0.15 


0.20 


Fig.  4.  Comparison  of  the  center  deflection  predicted  by  the  analytical  and  finite  element 
methods  for  the  case  of  symmetric  cross-ply  CFRP  laminates  (m,90,0,90,0)s  . 


(a) 


Time  (sec) 


(b) 


Time  (sec) 


Fig.  5.  Comparison  of  the  center  displacements  by  the  different  plate  theories  for  simply 
supported  cross-ply  CFRP  laminates;  (a)  laminate  (m,90,0,90,0)s,  (b)  laminate  (0,90,0, 90, m)s  . 
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Table  2  Selected  center  displacement  values  on  the  symmetric  cross-ply  CFRP  laminates  for  the 
different  laminate  theories. 


Center  Displacement  ( m  ) 

(m,90,0,90,0)s  laminate  (90,0,90,0,m)s  laminate 


CLPT 

FSDT 

TSDT 

CLPT 

FSDT 

TSDT 

0.0005 

3.93E-04 

3.92E-04 

3.89E-04 

3.94E-04 

4.00E-04 

3.89E-04 

0.0010 

5.21E-04 

5.06E-04 

4.87E-04 

4.87E-04 

4.92E-04 

4.57E-04 

0.0015 

3.32E-04 

2.96E-04 

2.57E-04 

2.11E-04 

2.09E-04 

1.51E-04 

0.0020 

-3.72E-05 

-8.25E-05 

-1.20E-04 

-2.21E-04 

-2.30E-04 

-2.74E-04 

0.0025 

-3.52E-04 

-3.75E-04 

-3.85E-04 

-4.82E-04 

-4.90E-04 

-4.72E-04 

0.0030 

-4.36E-04 

-4.11E-04 

-3.73E-04 

-3.78E-04 

-3.75E-04 

-2.83E-04 

0.0035 

-2.54E-04 

-1.86E-04 

-1.18E-04 

1.18E-05 

2.37E-05 

1.35E-04 

0.0040 

6.08E-05 

1.34E-04 

1.91E-04 

3.88E-04 

4.01E-04 

4.38E-04 

0.0045 

3.15E-04 

3.46E-04 

3.47E-04 

4.68E-04 

4.70E-04 

3.81E-04 

0.0050 

3.61E-04 

3.22E-04 

2.61E-04 

1.94E-04 

1.81E-04 

1.35E-05 

Time  (sec) 


Fig.  6.  The  effect  of  the  lamina  material  properties  on  the  damping  of  deflection  in  symmetric 
cross-ply  laminates  (m,90,0,90,0)s . 


The  deflection  suppression  time  is  defined  as  the  time  required  to  reduce  the  uncontrolled 
center  deflection  to  one-tenth  of  its  magnitude.  The  deflection  suppression  time  ratio 
(suppression  time  divided  by  the  maximum  suppression  time)  can  be  shown  to  be 
Ts  =  hm  /2 zm  ,  where  hm  is  the  thickness  of  the  magnetostrictive  layer  and  zm  is  the  positive 
distance  between  the  mid  plane  of  the  magnetostrictive  layer  and  the  mid  plane  of  the  plate. 
The  maximum  deflections  (JVmax)  of  the  composite  plate  and  the  suppression  times  for  the 

different  position  of  smart  layers  are  presented  in  Table  3.  The  effect  of  the  smart  layer 
positions  on  the  vibration  suppression  can  be  shown  in  the  Figure  7.  It  is  observed  that  as  the 
smart  material  layer  is  moved  farther  from  the  mid-plane  the  suppression  time  decreases,  as 
may  be  expected  because  of  the  moment  effect  by  smart  layer  actuations. 
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Time  (sec) 

Fig.  7.  The  effect  of  the  smart  material  layer  position  on  the  deflection  for  the  symmetric  cross- 
ply  CFRP  laminates. 

Table  3  Vibration  suppression  time  for  the  different  smart  layer  positions  on  the  symmetric 
cross-ply  CFRP  laminate  (m,90,0,90,0)s 


Lamination 

Scheme 

z,n  (m) 

Ts 

Wmzx  (10  4  W  ) 

t  at  ^ax/ 

t  at  /1Q 

(m,90,0,90,0)s 

0.045 

0.11 

1 

5.21 

0.0285 

(90, m, 90, 0,90) 
s 

0.035 

0.14 

3 

5.09 

0.0350 

(0,90,m,90,0)s 

0.025 

0.20 

0 

4.90 

0.0480 

(90, 0, 90, m, 90) 
s 

0.015 

0.33 

3 

4.85 

0.0850 

(0,90,0, 90,  m)s 

0.005 

1.00 

0 

4.87 

0.2560 

The  effect  of  the  thickness  of  smart-material  layer  on  deflection  damping  characteristics  is 
studied  next.  It  is  observed  that  thicker  smart  material  layers  result  in  better  attenuation  of 
the  deflection.  This  is  due  to  a  larger  mass  inertia  that  is  caused  by  the  large  increase  in  the 
moment  of  inertia  of  the  system  when  thickness  of  the  smart  material  layer  is  increased.  We 
note  that  the  smart  material  layer  has  a  density  of  five  times  that  of  the  composite  material. 
The  suppression  times  and  characteristics  for  different  smart  layer  thicknesses  are  shown  in 
Table  4  and  Figure  8. 

Figure  9  shows  the  effect  of  the  feedback  coefficient  c(t)kc  on  the  vibration  suppression 

characteristics.  Two  different  values  of  the  feedback  coefficient  are  used;  104  and  1 0  .  It  can 
be  seen  that  the  suppression  time  increases  when  the  value  of  the  feedback  coefficient 
decrease.  This  is  expected  because  the  coefficients  of  the  damping  matrix  decrease,  thereby 
resulting  in  less  damping. 
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Table  4  Suppression  times  for  the  different  smart  layer  thicknesses  in  symmetric  cross-ply 
laminates  (m,90,0,90,0)s 


Lamina  Thickness 
( mm ) 

Z,n  (fn) 

Ts 

W 

max 

( 1 0  4  m ) 

W  / 

t  at  ™x/o 

K 

=  io,  hm  =  2 

0.0410 

0.0244 

4.77 

0.0780 

K 

=  io,  K  =  4 

0.0420 

0.0476 

4.98 

0.0485 

K 

=  10,  hm  =  5 

0.0425 

0.0588 

5.05 

0.0420 

K 

=  10,  K  =  6 

0.0430 

0.0698 

5.10 

0.0350 

K 

=  10,  K  =  8 

0.0440 

0.0909 

5.17 

0.0320 

K ■ 

=  10,  K  =  10 

0.0450 

0.1111 

5.21 

0.0285 

K 

=  5,  hm  =  5 

0.0225 

0.1111 

9.12 

0.0310 

K 

=  5,  hm  =  10 

0.0250 

0.2 

9.45 

0.0400 

Time  (sec) 


(b) 


Time  (sec) 


Fig.  8.  The  effect  of  the  thickness  of  smart  material  layers  on  the  deflection  damping 
characteristics  of  symmetric  cross-ply  laminates  (m,90,0,90,0)s . 
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3 


-»-c(t)Kc=1000 


Fig.  9.  Effects  of  the  magnitude  of  the  feedback  coefficients  on  the  suppression  time  for 
symmetric  cross-ply  CFRP  laminates  (m,90,0,90,0)s . 

The  deflection  damping  characteristics  of  symmetric  angle-ply  and  general  angel-ply  composite 
laminates  are  studied  using  full  plate  F.E.  models.  Observations  made  earlier  on  various 
characteristics  such  as  the  effects  of  smart  layer  position,  its  thickness,  and  magnitude  of  the 
feedback  coefficient,  are  also  valid  for  these  laminates,  as  shown  in  Figures  10  and  1 1. 


.  Time  (sec) 

(b) 

Fig.  10.  Center  displacement  versus  time  for  symmetric  angle-ply  CFRP  laminates. 
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0.0006 


- (m,45,-45,0,90)s 


(a) 


Time  (sec) 


0.0004 


0.0002 


0.0000 


Q  -0.0002 


-0.0004 


0.0 


(b) 


(45,-45, 0,90, 30, -30, m,60)s 
(45,-45, 0,90, 30,  m, -30, 60)s 
-  (45,-45,0,90,m,30,-30,60)s 


o.i 


0.2  0.3 

Time  (sec) 


0.4 


0.5 


.  ,  Time  (sec) 

(c) 

Fig.  1 1 .  The  effect  of  the  position  of  the  smart  layer  on  the  deflection  of  symmetric  general 
angle-ply  CFRP  laminates. 


5.3  Clamped  composite  laminates 

Next,  fully  clamped  laminated  plates  are  analyzed  using  8x8  mesh  in  a  full  plate.  The  effect  of 
the  boundary  conditions  on  the  deflection  is  shown  in  the  Figure  12.  The  maximum 
displacements  of  the  simply  supported  plate  are  greater  than  those  of  the  clamped  case,  which 
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is  expected.  Simply  supported  laminates,  which  have  larger  displacements,  take  less 
suppression  time  compared  to  the  clamped  laminates. 


Time  (sec) 


Time  (sec) 

Fig.  12.  The  effect  of  boundary  conditions  on  the  center  displacements  for  the  different 
laminates;  (a)  Comparison  of  simply  supported  and  clamped  laminates,  (b)  results  for  different 
smart  layer  position. 

5.4  Effects  of  the  mechanical  loading 

Numerical  studies  are  also  carried  out  to  analyze  smart  composite  laminates  under  uniformly 
distributed  load  q0  instead  of  specified  initial  velocity  field.  Figure  13  shows  the  center 

deflection  for  selected  simply  supported  and  clamped  laminates  under  continuously  applied 
uniformly  distributed  loading,  while  Figure  14  shows  the  case  under  suddenly  applied  step 
loading.  The  effect  of  sinusoidal  loading  on  the  central  displacement  has  been  studied.  The 
results  of  symmetric  cross-ply  laminates  with  simply  supported  boundary  conditions  and 
subjected  to  sinusoidal  and  uniformly  distributed  loads  are  shown  in  Figure  15.  The 
transverse  displacements  in  Figures  13-15  are  plotted  in  the  nondimensionalized  forms  as 
F  h 3 

w0  a  oo) -i— . 

b  q  0 
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Fig.  13.  Nondimensionalized  center  deflection  versus  time  for  simply  supported  and  clamped 
laminated  plates  under  uniform  load. 


Time  (sec) 


Fig.  14.  Nondimensionalized  center  deflection  versus  time  for  simply  supported  and  clamped 
laminated  plates  under  suddenly  applied  uniform  load. 


Time  (sec) 


Fig.  15.  Comparison  of  center  displacement  versus  time  for  simply  supported  laminated  plates 
under  sinusoidal  and  uniformly  distributed  loads. 
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6.  NONLINEAR  FINITE  ELEMENT  MODEL 
6. 1  Virtual  Work  Statements 


The  weak  forms  of  the  nonlinear  equations  of  motion  are 


0  = 


-Nrr+—^N„+Su  o 


Jn‘  ^  dx  dy 

{^o (NxxrJx+Nxyr]y)}ds 


r  d2u0  1  r  d2(/>x  ~  ,  „ 

dr  dr  dt2 


2  ^3w,  ^ 
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xy 


dy 
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where  V  n  and  p  are  defined  as 
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d\Vn 

The  primary  variables  of  the  third-order  theory  are  un,  us,  w0,  — — ,  <j>n,  <j>s ,  where  ( un,us )  denotes 

dn 

in-plane  normal  and  tangential  displacements,  and  (<j>n,<j>s)  are  the  rotations  of  a  transverse  line 
about  the  in-plane  normal  and  tangent.  Lagrange  interpolation  functions  of  (un,us,<f)n,<f)s)  and 
Hermite  interpolation  function  for  w0  are  used  for  the  formulation  of  the  displacement  finite 

element  model.  A  conforming  element  that  has  eight  degrees  of  freedom 

(u0,v0,w0,  w0x ,  w0  ), ,  w0 xy , <j>x ,  <j)y )  is  used. 
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where  y/^  denotes  the  Lagrange  interpolation  functions  and  (pet  are  the  Hermite  interpolation 
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This  completes  the  nonlinear  finite  element  model  development  of  the  third-order  shear 
deformation  plate  theory.  As  in  the  case  of  linear  finite  element  analysis,  the  nonlinear 
transient  equations  may  be  solved  using  the  Newmark  time  integration  schemes. 


7  CONCLUSIONS 

A  unified  third-order  plate  theory  is  used  to  develop  a  displacement  finite  element  model  of 
laminated  composite  plates  with  smart  material  layers  for  transient  response.  The  unified 
formulation  includes  the  classical,  first-order,  and  third-order  shear  deformation  plate  theories 
as  special  cases.  The  smart  material  used  in  this  study  to  achieve  damping  of  transverse 
deflection  is  the  Terfenol-D  magnetostrictive  material,  although  in  principle  any  other  actuating 
material  can  be  used.  The  negative  velocity  feedback  control  is  used  to  damp  the  transient 
response.  A  number  of  parametric  studies  were  carried  out  to  understand  the  damping 
characteristics  of  various  laminates  with  embedded  smart-material  layers.  Some  of  the 
observations  are  summarized  below. 

Use  of  quarter  plate  models  in  place  of  full  plate  models  was  studied  first.  It  is  found  that  for 
antisymmetric  cross-ply,  angle-ply,  and  general  angle-ply  laminates  and  symmetric  cross-ply 
laminates  with  simply  supported  boundary  conditions,  a  quadrant  of  the  plate  with  proper 
symmetry  boundary  conditions  may  be  used  to  reduce  the  computational  effort. 

The  maximum  damping  of  deflection  occurs  when  the  smart  material  layers  are  placed  farthest 
from  the  midplane  because  the  bending  moments  generated  by  the  smart  material  layer  are 
maximum  when  they  are  away  from  the  midplane.  It  is  observed  that  for  a  lower  value  of  the 
feedback  coefficient  the  time  taken  to  damp  deflection  is  longer.  The  damping  tendencies  of 
symmetric  cross-ply,  angle-ply  and  general  angel-ply  composite  laminates  have  been  found  to 
be  similar.  The  damping  characteristics  of  a  fully  clamped  composite  plate  are  similar  to  that 
of  a  simply  supported  composite  plate. 
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